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Abstract 

Spatial resolution of less than 200 is challenging for thermal neutron detection. A novel readout scheme based on the time- 
projection-chamber (TPC) concept is used in a gaseous electron multiplier (GEM) detector ||T|. Thermal neutrons are captured in 
a single 2 jum thick Boron-10 converter cathode and secondary Helium and Lithium ions are produced with a combined energy 
of 2.8 MeV. These ions have sufficient energy to form straight tracks of several mm length. With a time resolving 2-dimensional 
readout of 400 fim pitch in both directions, based on APV25 chips, the ions are tracked and their respective origin in the cathode 
converter foil is reconstructed. Using an Ar-C02 93:7% gas mixture, a resolution of 100 jum (FWHM 235 //m) has been observed 
with a triple GEM-detector setup at the Garching neutron source (FRMII) for neutrons of 4.7 A. 
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1. Introduction 

Thermal neutron detection with spatial resolution around or 
below 0.2 mm will become mandatory for modern thermal neu¬ 
tron crystallography or fast imaging detectors 0. 

Gaseous micro-pattern detectors like GEM (gaseous electron 
multiplier) or Micromegas have proven to reach spatial resolu¬ 
tions well below 100 jim 0, depending on the fineness of their 
segmentation. Here a self calibrating TPC-like analysis for this 
kind of detector is presented with a GEM-detector. 

2. TPC-like Analysis of Ion Tracks 

A triple GEM-detector with segmented 2D-readout of 0.4 
mm pitch was used in combination with a ^°B converter cath¬ 
ode to detect thermal neutrons. Thermal neutron capture in '‘'B 
has a high cross section and produces a ^Li and a ^He ion back- 
to-back with a combined energy of 2.8 MeV. 

-H n Li He + y+ Q{2.19MeV){9A%) 

-H n Li He + Q{2.19MeV){6%) 

The heavy He or Li ion create a straight line of ionized elec¬ 
trons in the drift-volume of a micro-pattern detector. The full 
track information can be reconstructed from the electron drift 
time measured as a function of the responding readout seg¬ 
ments. This allows for reconstruction of the neutron interaction 
point in the converter cathode. 

The APV front-end-electronics samples the signal on every 
strip in steps of 25 ns. With typical rise times of 150 ns the 
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beginning of a signal, and thus the drift time of the respective 
electrons, can be determined with a fit to the rise of the charge 
distributions. 



Figure 1: The starting time of the signal is determined for every strip by fit to 
the rise time of the respective charge signal. Direction and orientation of the 
track are determined by a linear fit to the start points 

From a fit to the starting times as a function of hit strips the 
inclination angle of the track can be reconstructed. 

3. Determination of spatial resolution 

As conversion cathode a 2 fim layer of ^°B evaporated on an 
aluminum supporting sheet has been used. It was mounted to a 
triple GEM detector with two-dimensional strip readout, based 
on crossed strips with a pitch of 400 /urn in both directions. The 
readout anode was delivered by CERN. X and Y directions are 
covered by 256 strips each with a copper strip width of 80 fim 
for the upper layer and 320 ijm for the lower layer. All the insu¬ 
lator material but the direct insulation between X and Y copper 
strips have been etched off, such that the different copper width 
guarantees similar pulse-height in both planes. The readout is 
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Figure 2: Experimental setup for position resolution measurements at the 
Garching Research Reactor (FRMII) at the beam line TREFF. Neutrons of 4.7A 
(13.4 meV) are collimated through a system of four apertures to a slit of 250 
/im width and 65 mm height on the converter cathode of a GEM-detector 


based on APV-25 front-end boards |j4l and the SRS data acqui¬ 
sition system Q and was triggered by the signal on the lower 
side of the last GEM-foil, which is in principle the inverted sig¬ 
nal of the readout plane. The drift gas, Ar-C02 of 93:7 Vol.% 
at atmospheric pressure, dehnes in combination with the thick¬ 
ness of the absorber the maximum range of the ^He ions to 12 
mm. The width of the drift gap of 6 mm was chosen to decrease 
the impact of the Bragg-peak in the track reconstruction and 
to assure that only tracks with large inclination angle stop in 
the active volume. The distances between GEM-foils and the 
distance between last GEM-foil and the readout anode was 2 
mm each. Eor optimized spatial resolution the drift velocity of 
the electrons was chosen quite low to 11 jum/ns, by adapting 
the electric held in the drift region to the time resolution of the 
APV25 chips. 

Neutrons of 4.7 A are available at the Garching ERMll re¬ 
search reactor at the TREEE beam line. The beam was colli¬ 
mated by a system of four apertures, see hgure|^ The collima- 
tion was done by the collimators 1, 2 and 4 and a possible beam 
halo was cut by collimator 3. The beam was formed to a slit of 
180 jum width and 65 mm height, dehned by the collimators on 
the last aperture, which lead to a maximum width of the beam 
spot on the converter cathode of 250 fim with only 0.32 mrad 
divergence. To determine the neutron prohle on the conversion 
layer the hit distribution has been approximated by a boxcar 
function. The convolution with the Gaussian-distribution de¬ 
scribing the spatial detector resolution allows then to simulate 
the measured distribution. This leads to a double error-function 
of the form: 


fix) = A er/ 




+ erf 


-jx-fJ.) 
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cr 
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A denotes a normalization factor, a is the half width of the 
slit, jj. is the center of the slit and cr is the standard deviation of 
the Gaussian, which represents the detector-resolution. 

In a standard analysis the neutron position would be deter¬ 
mined by the center of the charge distribution for every track in 
the readout plane. The error in position resolution is therefore 
limited by the inclination angle and the length of the track of 
the secondary ion, see figure The resolution achievable by 
this method is a EWHM of 3.4 + 0.1 mm. 

The resolution can be improved by determining the starting 
point of the secondary ion track and thus the point of conversion 
in the cathode. 




Ax 


Figure 3: The conversion point in the cathode can be determined starting from 
the center of charge position, which is well defined in the track, and applying a 
con'ection to the position, which is dependent on the direction and length of the 
track and the type of the ion producing the track 



Figure 4: Projection onto the X-direction for charge averaged centers of tracks 
of a slit for a slit of 180 /im width. The width of the distribution is 3.4 mm 
(FWHM) (o- = 3.667x0.4 mm ). 
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Figure 5: Projected track length in X-readout direction plotted against the cen¬ 
ter of charge position for particles with negative inclination angle going from 
right to left in X-direction. A linear fit is applied to compensate for the track 
length dependence. 

With a single GEM-detector spatial resolution of below 
470 yum could already be achieved by determination of the last 
strip, which was hit in a track(||6l). In a triple GEM-setup the 
spread of the charge cloud by diffusion towards the anode is 
quite significant. This can be overcome by fully tracking the 
ions in the drift volume. 

The method would work best if all ions would transverse the 
complete drift region and if the correlation between responding 
readout strips and ion track would be precise. Certainly, ions 
with large inclination angles stop within the 6 mm wide drift re¬ 
gion. The systematic effects: electron diffusion and capacitive 
coupling of neighboring readout strips, lead to an overestima¬ 
tion of the length of the track on both ends for both traversing 
and stopping ions. Therefore, the TPC-like reconstruction is de¬ 
pendent on the track length, the inclination angle, and the non- 
homogeneous the charge distribution in the track, the skewness. 
Skewness and track length are dependent on the particle type as 
well. 

In the following the dependencies of the reconstruction on 
the track length, and the skewness of the charge distribution in 
the track will be discussed. These corrections are performed by 
linear fits to internal detector correlations see figure]^ 

Eigure shows a clear correlation between the projected 
track length in one readout direction of a track and the recon¬ 
structed position. The length of the track is determined from a 
Gaussian fit to the charge distribution on all responding strips 
individually in X and Y direction. The dependence of the cen¬ 
troid position on the projected track length was determined by 
linear fit and the position was corrected depending on the pro¬ 
jected track length. 

After application of this correction the dependence of the re¬ 
constructed position on the skewness of the charge distribution 



Figure 6: Relative position of the center of charge in the track plotted against 
the track length corrected position. A center of charge position at the end of the 
track corresponds to 1 and (-1) to the begining of the track. The dependence is 
determined by a linear fit and the position is corrected for this dependence. 
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Figure 7: Projected angle in X-readout direction plotted against corrected posi¬ 
tion. A linear fit is applied to compensate for the angular dependence. 
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Figure 8: Hit pattern in one readout direction, after applying correction factors 
on track length, angle and charge distribution in the track. The achieved res¬ 
olution is determined by a Gaussian function convoluted with a double error 
function and was determined to (235± 20) //m FWHM (cr = 0.252x0.4 mm ) 

is shown in figure The linear fit gives the parameters for 
the correction in a second step and in the last step the fit to the 
reconstructed angle of the track as a function of reconstructed 
position is investigated, as shown in figure]^ 

Since there are dependencies of all these measurement pa¬ 
rameters it is not advantageous to correct them iteratively in a 
calibration one after the other. 

The result of this calibration is shown in figure where the 
projection of the strip in the X-direction is shown and the reso¬ 
lution is determined by a fit of equation A spatial resolution 
of (235 + 25) fim FWHM (cr =(100 + 15) yum) is achieved. 

4. Conclusion 

Using juTPC analysis the inclination angle of the track is de¬ 
termined, which defines the direction of the track. As origin the 
centroid of the tracks was determined and its dependence on the 
track length, inclination angle, and the skewness of the charge 
distribution in the track is examined. The iteration of calibration 
factors determined from linear fits to the reconstructed position 
in dependence to the mentioned track parameters were adapted 
to compensate for these. Altogether this corrections in addition 
to a cut on the minimal track length of 3.5 strips in both direc¬ 
tions led to a gaussian distribution with an improved resolution 
of (235+25) jdm FWHM (cr = 100+15 yum) in comparison to a 
standard centroid determination with a resolution of 3.4 mm. 

This procedure produces implicitly correction factors that 
have been so far obtained in MIP tracking experiments by ded¬ 
icated calibration runs with known incidence angle of the parti¬ 
cle beam. 
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